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irradiation: Interplay of geometry and dewetting 



S. Afkhami and L. Kondic 

Department of Mathematical Sciences and Center for Applied Mathematics and Statistics, 
New Jersey Institute of Technology, Newark, NJ 07102 USA 

Metallic nanoparticles, liquefied by fast laser irradiation, go through a rapid change of shape 
attempting to minimize their surface energy. The resulting nanodrops may be ejected from the 
substrate when the mechanisms leading to dewetting are sufficiently strong, as in the experiments 
involving gold nanoparticles [Habenicht et al.. Science 309, 2043 (2005)]. We use a direct continuum- 
level approach to accurately model the process of liquid nanodrop formation and the subsequent 
ejection from the substrate. Our computations show a significant role of inertial effects and an 
elaborate interplay of initial geometry and wetting properties: e.g., we can control the direction 
of ejection by prescribing appropriate initial shape and/or wetting properties. The basic insight 
regarding ejection itself can be reached by considering a simple effective model based on an energy 
balance. We validate our computations by comparing directly with the experiments specified above 
involving the length scales measured in hundreds of nanometers, and with molecular dynamics 
simulations on much shorter scales measured in tens of atomic diameters, as in M. Fuentes- Cabrera 
et al., Phys. Rev. E 83, 041603 (2011). The quantitative agreement, in addition to illustrating how 
to controlling particle ejection, shows utility of continuum-based simulation in describing dynamics 
on nanoscale quantitatively, even in a complex setting as considered here. 

PACS numbers: 47.11.-j,68.08.-p,81.16.Rf,68.55.-a 



Evolution of fluid drops deposited on solid substrates 
has been a focus of large research effort for decades. More 
recently, this effort has been particularly extensive on 
the nanoscale, due to relevance of nanostructures in a 
variety of fields, ranging from DNA sequencing to plas- 
monics [l|, [H , nanogmagnetism fl| , and liquid crystal dis- 
plays and solar panel designs [4]. For example, the size 
and distribution of metallic particles strongly affects the 
coupling of surface plasmons to incident energy [1] . Con- 
trolling this coupling has the potential for large increases 
in the yield of solar cell devices. 

In addition to physical experiments, modeling and sim- 
ulations provide a significant insight into the effects gov- 
erning the evolution of drops and other liquid structures. 
Numerical simulation has an advantage of allowing to 
switch off some of the involved physical phenomena and 
therefore isolate the dominant ones. On nanoscale, how- 
ever, it is not trivial to decide on appropriate simula- 
tion technique. Molecular dynamics (MD) simulations, 
while very powerful and presumably as close to real phys- 
ical picture as it is our knowledge of underlying interac- 
tion laws, are still extremely computationally demanding. 
Therefore, one would like to resort to the continuum- 
based simulation. However, it is not immediately clear 
that this approach is appropriate on the nanoscale where 
the basic assumptions of continuum fluid mechanics are 
pushed to their limits. In addition, the continuum simu- 
lation of free surface time-dependent problems based on 
full three-dimensional Navier- Stokes (N-S) formulation 
is still computationally demanding. For this reason, this 
type of simulation is rarely attempted in practice and 
instead the researchers resort to asymptotic (long-wave) 
methods. Such an approach has been applied with suc- 



cess to problems involving wetting and dewetting of drops 
and films; see ^ for reviews. This direction is, however, 
questionable due to the inherent assumptions, in partic- 
ular when contact angles are large and inertial effects are 
significant. One such class of problems involves liquid 
metals. These configurations may even lead to topolog- 
ical changes under liquefaction with metal particles de- 
taching from the substrate [6]. In such cases, long- wave 
theory clearly cannot be applied. 

In this work, we show that direct numerical solutions 
of N-S equations lead to the results which are in quan- 
titative agreement with (i) MD simulations of dewetting 
of 1 — 1.5 nm thick, liquid copper (Cu) disks of radius of 
10—15 nm, and (ii) physical experiments involving 'jump- 
ing' gold (Au) nanoparticles of typical in-plane length 
scale in the range of 100 nanometers. While utility of 
continuum mechanics on nanoscale problems involving 
liquid metals was discussed previously (see, e.g., [3]), as 
far as the authors are aware, this is the first attempt to 
compare explicitly continuum N-S simulations both to 
MD and to physical experiments in a dynamic setting. 
Furthermore, we utilize the simulations to illustrate the 
dominant effects and discuss in particular the conditions 
leading to ejection of fluid material from the substrate. 
We also present an effective model based on an energy 
balance to provide a basic insight into the ejection mech- 
anism. We first consider dewetting of liquid Cu disks 
recently simulated by MD Q and then shift focus to Au 
particles considered experimentally [g] . 

Methods. The simulations are based on the Volume 
of Fluid (VoF) approach, which, when coupled to a flow 
solver, is also used to compute quantities related to sur- 
face tension that enter the flow calculation. Until re- 
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cently, the VoF method was not deemed appropriate for 
the study of phenomena in which surface tension is the 
driving force. Recent improvements to calculating cur- 
vature and applying the surface tension force appear to 
resolve this issue [9|. Within this approach, one solves 
the three-dimensional N-S equations that govern the mo- 
tion of the fluid both inside of the liquid domain, and 
in the surrounding gas phase. For Cu disk structures, a 
free-slip boundary condition is specified at the substrate; 
this choice is motivated by the MD simulations |8j. This 
assumption is also consistent with the fact that the length 
scale associated with slip is nanometric for most sys- 
tems lip|. When considering Au structures, the Navier 
slip [Uj with slip length of 3 nm is imposed; this par- 
ticular value is obtained by direct comparison with the 
experiments Q and is further discussed in Supplemen- 
tary Materials [l2|. Such a value is also consistent with 
the MD simulations [l3|, indicating slip lengths of a few 
nanometers for mesoscopic systems. The physical prob- 
lem considered is the one of a Newtonian, isothermal, in- 
compressible fluid and therefore in the present work we 
do not consider the effects related to phase change and/or 
thermal variation of material properties occurring in ex- 
periments Q . Simulations are carried out on an adaptive 
mesh [l^, and numerical convergence has been verified 
by considering simulations with different grid resolutions 
See [ill for more details. 

Results. Figure [TJa) shows the evolution of a Cu disk 
of initial height ho = 15 A and radius = 150 A with 
an initial contact angle Oq = 90°, when the equilibrium 
contact angle ^eq = 80°, that dewetts and collapses into 
a spherical cap. Figure [TJb) and (c) show the front posi- 
tion, R{t)^ and velocity, R{t)^ respectively. We also show 
the results of the MD simulations obtained by employ- 
ing two different Lennard- Jones (LJ) potentials, LJ-a and 
LJ-b, that differ by the depth of potential well [8|. These 
results show that our numerical results and the MD sim- 
ulations yield fully consistent retraction time scales. We 
also see that 6>eq influences the initial rapid retraction, 
with larger 6>eq leading to faster dynamics for early times. 

A different type of evolution occurs if ^eq is further in- 
creased. To illustrate this. Fig. [2fa) shows snapshots of 
the profiles for ho = 15 A, Rq = 150 A, Oq = 90°, and 
^eq = 140°. For this ^eq^ the contraction is so fast that 
the nanodrop jumps off the surface, following elongation 
in the y direction (normal to the substrate). Similar be- 
havior is observed for any 6>eq ^ 130°, again in quantita- 
tive agreement with the MD simulations [8|. 

To gain additional insight regarding the effects that 
drive contraction, and possibly ejection of the fluid, we 
have carried out additional set of simulations with mod- 
ified inert ial effects. This was implemented by vary- 
ing incompressible fluid density, but more generally the 
dynamics can be characterized by Ohnesorge number. 
Oh = r]/^/ paRo , which is related to the Reynolds number 
defined based on a characteristic capillary velocity, a/r], 
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FIG. 1: (Color online) Evolution of a Cu disk {ho = 15 A, 
Ro = 150 A, Oo = 90°); (a) Snapshots for 6'eq = 80°; (b) 
R{t) and (c) R{t) for (■) 6>eq = 80° and (a) 6'eq = 115°. MD 
simulations for (•) LJ-a, Oeq = 116°, and (♦) LJ-b, ^eq = 
75.6°, with ho = 10 A and (a) LJ-b with /lo = 15 A 0. 
Predictions of the model, Eq. ([T]), for ^eq = 80° (solid green 
line) and 115° (dashed blue line). 

as Re = Oh~'^. For liquid Cu, Oh ^ 0.35, suggesting 
that inertial effects are important; the same conclusion 
can be reached by considering an intrinsic length scale, 
£v = rf /{p(j)^ above which inertial effects become signifi- 
cant; for liquid Cu, ^ 1.78 nm, therefore smaller than 
typical length scales considered here. As an example. 
Fig. [2];b) and (c) shows R(t) and R(t) for Oh ^ 0.35 and 
Oh ~ 1.1. Reduced inertial effects eliminate the most 
noticeable feature of the dewetting: nanodroplets do not 
detach from the surface for the same ^eq- Additional sim- 
ulations have shown that even for ^eq = 150° there is no 
detachment for Oh ^1.1. 

A basic understanding of the effects that drive the fluid 
evolution can be reached by a relatively simple effective 
model based on the balance of relevant energies. In this 
model, the evolution of the nanodisk obeys the energy 
balance equation 



(1) 



where Ek, Eg are the kinetic and surface energy, re- 
spectively, and D is the rate of energy loss due to 
viscous dissipation, neglecting the gravitational energy. 
In this model, we describe the nanodisk dynamics by 
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FIG. 2: (Color online) Evolution of a disk {ho = 15 A, 
Ro = 150 A, 6*0 = 90°); (a) for 6>eq = 140°, the disk detaches 
at around 100 ps (see also the animation in p^); (b) R{t) 
and (c) R{t) for 6'eq = 80° (■)/( ) and 140° (•)/( ) when 
Oh ^ 0.35 (black) and Oh ^1.1 (green). Predictions of the 
model, Eq. Q, for 6'eq = 80° (solid red line) and 6'eq = 140° 
(dashed blue line) when Oh ^1.1. 



a fluid cylinder evolving on a solid substrate; a simi- 
lar model has been considered for a drop impact prob- 
lem [l^. Here, = p J^\V\'^dQ/2 where the integra- 
tion is over the fluid cylinder, and V = V{x,t) is the 
axially symmetric fluid velocity. Using Young's law [Tsj , 
Es=cr [7TR{t)^{l - cos6>eq) + 27TR{t)h{t)] , whcrc R{t) is 
the wetted radius and h{t) is the cylinder height. The 
viscous dissipation energy D = J^if : \/V)d^}, where 
f is the shear stress tensor. Using axisymmetric stag- 
nation point flow and the conservation of mass, Eq. ([T]) 
reduces to a nonlinear second-order variable-coefficient 



ordinary differential equation given in [12|. This equa- 
tion is then solved numerically, and the results are shown 
together with the computational results in Figs. [1] and [2] 
(plane lines). We find very reasonable degree of agree- 
ment, suggesting that this simple model captures well the 
main mechanisms driving fluid contraction. 

Figure [3] shows additional comparison between the nu- 
merical simulation and the effective model. In Fig. [3l^a), 
we observe that both the simulations and the model, 
Eq. ([T]), predict approximately linear dependence of the 
ejection velocity on (1 — cos6>eq). We expect that the dis- 
crepancy between the two is maily due to the deficiency 
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FIG. 3: (Color onhne) (a) The effect of ^eq on the ejection 
velocity, Vd, for a Cu nanodisk {ho — 15 A, Ro = 150 A) 
predicted by the model based on Eq. {[T]) (a), see also Eq. (13) 
in |T2], and simulations (■). (b) Phase diagram showing the 
influence of Oh and ^eq on the ejection. The solid line is 
based on the model, Eq. ([1]), see also Eq. (13) in |lj], and the 
symbols mark the result of the numerical simulation showing 
detaching (■) and non-detaching (a) nanodrops. 
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FIG. 4: Evolution of a Au equilateral triangle (a = 405 
nm, ho = 47 nm, 6*0 = 90°, and 6'eq = 140°). The structure 
collapses into a droplet and detaches at t ^ 7 ns (see also the 
animation in [l3lV 



of the model for large 6'eq [l2| . Figure ^h) show a phase 
diagram in (O/i, 6'eq) parameter space ihustrating the 
criteria for ejection resulting from the model, together 
with the results of the numerical simulation, concentrat- 
ing here on small values of Oh. Although the model 
underpredicts ^eq required for ejection, we note a consis- 
tent trend of the results. Note that for large Oh ^1.1 
considered previously, both the model and the simulation 
are in non-detaching regime. 

Next we consider the configuration representative of 
physical experiments where Au triangular structures 
were liquefied and let to evolve on Si02 substrate 
Figure |4] shows the snapshots of initially equilateral tri- 
angle (^0 = 90°). The dewetting process first starts at 
the vertices, where the curvature is high. Due to high 
surface tension and large ^eq, the fluid starts to accu- 
mulate there. The humps at the vertices then coalesce 
into a droplet. Owing to a low viscosity of liquid gold. 
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FIG. 5: Evolution of a Au isosceles triangle for ho = 24 nm, 
Oo = 90°, long sides a = 438 nm and short side b = 247 nm; 
^eq = 140° . The triangle collapses into a droplet and detaches 
from the substrate at around 5 ns (see also the animation 
in [13). 

inertial effects dominate over viscous dissipation (here 
Oh ^ 0.047 [lij), giving rise to an upward movement 
that leads to droplet detaching from the surface with a 
velocity of ~ 24 m s~^. This process is consistent with 
the dewetting induced ejection mechanism outlined in [6|. 
We note the consistency of the time scales found by the 
simulations and the experiments in [6[, where an ejec- 
tion time scale of the order of 10 ns was observed. In 
addition to comparing favorably with experiments, the 
simulations allow for an additional insight since it is pos- 
sible to explore very fast time scales, and also easily vary 
the parameters, such as 6>eq and the initial shape, and ask 
what is their influence on the outcome. 

To illustrate the influence of geometry, we consider 
isosceles triangles. Figure [5] shows a perhaps unexpected 
result, that despite the fact that the vertex at the small- 
est angle collapses the fastest due to a higher curvature, 
the vertices at the larger angles arrive sooner to the center 
(future work should analyze how general this result is). 
This mismatch excites oscillations of the droplet trans- 
lating into a tumbling movement after the ejection. 

Next, we ask whether the ejection angle, a, is always 
7r/2. We find that a can be influenced for asymmet- 
ric collapse, and present two means for controlling it by 
modifying: (i) ^eq and (ii) the initial geometry. Figure 
[6ja-b) shows that increasing ^eq alters the no-ejection 
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FIG. 6: x-y cross sections of Au isosceles triangles for ho 
nm and — 90°; (a-b) a — 438 nm and b — 247 nm for (a) 
^eq = 130°, (b) 140°, showing the increase of ejection angle, a, 
as ^eq is increased (note that in (a), a = 0); (c-d) Oeq = 145°, 
and (c) a = 362 nm, b = 185 nm, (d) a = 475 nm, b = 308 
nm; t = 0, 0.6, 1.2, 3, 4.4, 11, 19, 35 ns. (e-f) Ejection 
angle a versus the Oeq (e), when a/b ^ 1.77, and the initial 
geometrical ratio a/b (f), when ^eq = 145°. 



{a = 0) to a directional ejection (0 < a < 7r/2). A qual- 
itative understanding can be reached by recalling that 
a 7^ 7r/2 is due to the lack of synhronicity of the collapse. 
A larger ^eq results in higher surface energy, speeding up 
the fluid at the vertices. Therefore for larger ^eq, col- 
lapse is more synchronous and a is closer to 7r/2. For 
smaller ^eq, <^ may vanish (see Figure [6fa)), i.e. ejection 
does not happen, while an equilateral triangle with the 
same ^eq does eject, illustrating the effect of asymmetry. 
We next consider the aspect ratio of the lengths, (a, 6), 
of triangle sides, and find that the larger a/b ratios lead 
to increased asymmetry of the collapse and smaller a, as 
shown in Fig. [6]^ c-d). 

Fi gures ini^e) and (f) show o. versus ^eq for fixed a/b^ 
and versus a/b for fixed ^eq, respectively. In Fig. [6je) we 
see that the ejection angle increases monotonically with 
an increase of 6eq with no ejection for 6>eq ^ 135°. Figure 
[61(f) shows that a decreases with increased asymmetry 
and that for the considered 6>eq, the ejection does not oc- 
cur, i.e. a = 0, when a/b > 2.4. Very large ratios lead 
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to a diflPerent type of dynamics: Ray leigh- Plateau type 
of breakup, see [16| for discussion in the context of liquid 
metals. More elaborate analysis of the influence of geom- 
etry and wetting properties is needed to fully understand 
the dynamics - we leave this for future work. 

Conclusions. In this work, we have demonstrated that 
continuum simulations provide a good qualitative agree- 
ment with both MD simulations on the length scales in 
the range of 1 — 10 nm, and with the physical experi- 
ments with in-plane length scales measured in the range 
of 100 nanometers. We expect that this flnding will fur- 
ther motivate modeling and computational work on these 
scales, since it is suggesting that the continuum-based 
simulation has a predictive power. For the problem of 
dewetting and possibly detaching nanodrops, the simu- 
lations provide precise insight regarding the influence of 
inertial, viscous, and capillary forces, in addition to the 
liquid/solid interaction. This insight is also conflrmed 
by a simple model based on an energy balance and ac- 
counting for viscous energy losses. Furthermore, the sim- 
ulations provide a clear prediction that the direction of 
ejection of fluid from the substrate can be influenced (and 
controlled) by modifying either the wetting properties or 
the initial geometry. Future computational and modeling 
work will include the thermal and phase change effects, 
allowing to even more accurately model the dynamics of 
liquid metals on nanoscale. 
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SUPPLEMENTARY MATERIAL: 



COMPUTATIONAL MODEL 



MATHEMATICAL MODEL 

The equations of conservation of mass, \/ - U = 0, and 
momentum are written as 
d 



dt 



(pU) + V • {pUU) = -Vp -f^Fs, 



(2) 



where U = (ii, w) represents the velocity vector, p the 
pressure, p the fluid density, f the shear stress tensor, 
and Fst the surface tension force (per unit volume) acting 
on the fluid. The shear stress tensor is defined as f = 
piVU + iVU)^)^ where p represents the fluid dynamic 
viscosity. 

The flow equations have been written in an Eulerian 
frame of reference, and thus a solution of these equations 
must be coupled with a methodology for following the de- 
forming fluid-fluid interface. Here, the 'Volume of Fluid' 
(VoF) algorithm is implemented [l|-l3|. Volume tracking 
requires the introduction of a scalar function / defined 
as 

in fluid 1 
in fluid 2, 



/ 



for a two fluid system. Since / is passively advected with 
the flow, it satisfies the advection equation 



V)/ = 0. 



(3) 



Density and viscosity are then evaluated via volume- 
weighted formulae as p = pi + {p2 — pi)f and p = 
Pi + {p2 — Pi)fi respectively, where subscripts 1 and 2 
refer to fluids 1 and 2, respectively. 

The surface tension force in (|2]) is reformulated as 
an equivalent volume force [4], which is non-zero only at 
each interface cell 



(4) 



where a is a constant interfacial tension and 5s denotes 
the Dirac delta function for the surface separating the 
fluids. The curvature n and unit normal n directed into 
fluid 1 are geometric characteristics of the surface and 
are described in terms of / and computed with a second- 
order 'height-function' (HP) method Within the 
VoF-based sharp surface tension representation, dsn is 
equivalent to V/, thus 



Fst = (jnVf. 



(5) 



If partial slip is allowed at the contact line, we impose 
the Navier slip boundary condition [9] at the (static) solid 
surface y = 0, 



{u,w)\y=Q = Xd{u,w)/dy\y^Q, 



(6) 



where A is the slip length. For dynamic contact lines, the 
Navier slip boundary condition leads to a regularization 
of the viscous stress singularity at the contact line |9|]. 



The main issue is the imposition of the contact an- 
gle, described in [l, 0- Briefly, within our VoF frame- 
work, the contact angle boundary condition enters the 
N-S solver in two ways: it defines the orientation h of 
the VoF interface reconstruction at the contact line, and 
it influences the calculation of the contact line curvature; 
together, these two reflect the wettability of the surface. 
We implement the HF methodology to impose the con- 
tact angle in 3D The HF methodology accurately 
computes interface curvatures at the contact line, values 
that converge with mesh refinement. In our HF method- 
ology, the curvature at the contact line is computed at 
the desired contact angle in order to balance the surface 
tension forces on the contact line as in Young's relation 

Ii3 



Cr COS 6>eq = CFsg - CF si , 



(7) 



where Gsg and cFsi are interfacial tensions for the solid-gas 
and solid-liquid surfaces, respectively. If the actual con- 
tact angle is out of equilibrium, the resulting mismatch 
in curvature generates a kink in force thus driving the 
contact line to configuration that satisfies Young's rela- 
tion. 



EFFECTIVE MODEL 




FIG. 7: Cylindrical model of nanodisk retraction. V — 
[Vx^Vy) is the axially symmetric fluid velocity; x is the ra- 
dial coordinate and y is the height. 



Following [Uj, ll2| , we developed a simple model for the 
evolution of the nanodisk based on the energy balance. 
In this model, an ordinary differential equation for the 
wetted radius R{t) of a nanodisk is developed by model- 
ing the retraction dynamics of a fluid cylinder of radius R 
and thickness h evolving on a solid substrate, see Fig. [71 
In this model, the kinetic energy of a drop of volume 
is defined as 



Ek 



(8) 



The surface energy is defined as 



Es=(J [lTR(tf(l - cos6>eq) + 2lTR(t)h(t)\ . (9) 
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The rate of energy loss due to viscous dissipation is de- 
fined as 



D = 2ii 



j 

JQ. 



dx 



X 

1 fdV^^dV^ 

2 \ dy dx 



dy 



dn. 



We assume an axisymmetric stagnation point flow 



-2 dR{t) 
~R dt 



1 dR{t) 
R dt 



y 



1 dh{t) 
" hit)~dt 

-1 dh{t) 
2h{t) ~~dt 



-X. 



(10) 



(11) 



(12) 



We note that the velocity field assumed by Eqs. (pTj) and 
([12]) is shear free and so imply a slip at the solid sub- 
strate. This is consistent with the assumption of a free- 
slip boundary condition specified at the substrate for our 
nanodisk direct numerical solution. The energy balance, 
in dimensional form, reads 



d_ 

m 



(13) 



where we neglect the gravitational energy. Using the ve- 
locity field, Eqs. ([n]) and (O in Eqs. ([5]) and 1^0^, the 
volume conservation R{t)'^h{t) = RqHq^ and taking the 
time derivative in Eq. (fT3]) . we arrive at the following 
nonlinear second-order variable-coefficient ordinary dif- 
ferential equation 



A{R)R + B{R)R^ + C{R)R + D{R) = , 



(14) 



where 



A{R) 



'^^ 3 i?6 



B{R) 



4R^ 
B7 



D{R) 



2R 
Rih~o 



(1 -COS(9eq) - 



with the initial conditions R{t = 0) = Rq and R{t = 
0) = 0. We compute R{t) and R{t) by solving the initial 
value problem (fH|) using the forth-order Runge-Kutta 
method. One would naturally consider detachment when 
R = 0. Eq. (p!^ however becomes singular when R 
0. We propose the following procedure to alleviate this 
numerical difficulty. When detachment occurs, as t ^ 
oo, the total surface energy will be equal to the surface 
energy of a spherical droplet of radius Rs = (f^o^o)^^^ 
(see Fig. [7]), thus 



7TR{t oof {I - COs6>eq) + 27r 



R{t 



47rcr ( ^Rlho 



oo) 

2/3 



(15) 



Cu 



Au 



Ro (Cu), Lo(Au) 
A 

/)air(kg m~^) 

?7air(kg m"'' 
Pdrop(kg m~ 
?7drop(kg m"^s" 
a(kg s-2) 
Oh = Tj/yJ pcrL 



s 

-3\ 



150 A 
375 A 
2.93 A 
1.225 
0.00002 
7900 

0.004288 ((( 

1.305 

0.35 



405 nm 
494 nm 
1.5 nm 
1.225 
0.00002 
17310 
1500 K) 0.00425 {i 
1.15 
0.047 



1500 K) 



TABLE I: Overview of the sets of parameters used in the 
numerical simulations [Tsl . [T3 | . Initial conditions are specified 
in terms of disk radius, Rq for Cu and a typical length of a 
triangle side, Lo for Au. In computing Ohnesorge number, 
O/i, for the typical length scale, L, we use Rq for Cu and Lo 
for Au. 



where R{t oo) is the radius of the wetted base as 
t ^ oo. We use Eq. ([15]) to calculate R{t oo) as a 
function of 6>eq- Let us denote this detachment radius Rd- 
When numerically computed R{t) equals Rd^ we consider 
the drop detached and calculate its ejection velocity. If 
R{t) never reaches Rd^ we then consider the disk only 
collapsing into a sessile drop on the substrate. We note 
that Eq. (p!5]) leads to a real positive solution for Rd only 
for 6>eq < 107°, where R < 0.44. We use this value of R 
as the one at which the drop detaches from the substrate 
when (9eq > 107°. 



COMPUTATIONAL PARAMETERS 



Table |T] provid es an overview of the parameter sets 
used. The computational domain IS X Ly X Lr. and the 
maximum grid resolution of the adaptive mesh is repre- 
sented by A. For all the simulations, an open boundary 
condition (pressure and velocity gradient equal zero) is 
imposed at the top and a symmetry boundary condition 
is imposed on the lateral sides. 



DETERMINATION OF THE APPROPRIATE 
SLIP LENGTH 

We explored the influence of the slip length by carrying 
out numerical experiments for Au structures. Larger slip 
lengths lead to faster collapse and large detachment ve- 
locity. In Fig. m we show that by choosing the slip length 
of 3 nm, the computed detachment velocity of nanostruc- 
tures compares well with the experimental results from 
[l^ : this value is therefore used in all presented Au sim- 
ulations. 
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FIG. 8: Droplet detachment velocity for liquefied triangular 
Au nanostructures with a side length of 405 nm and vari- 
ous thicknesses, numerical results (■) and experimental mea- 
surements (■). Numerical results are obtained using the slip 
length of 3 nm. 
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